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Abstract 

At very high energy and wide angles, Compton scattering on the proton (7^ 
7p) is described by perturbative QCD. The perturbative QCD calculation has 
been performed several times previously, at leading twist and at leading order in 
as, with mutually inconsistent results, even when the same light-cone distribution 
amplitudes have been employed. We have recalculated the helicity amplitudes for 
this process, using contour deformations to evaluate the singular integrals over 
the light-cone momentum fractions. We do not obtain complete agreement with 
any previous result. Our results are closest to those of the most recent previous 
computation, differing significantly for just one of the three independent helicity 
amphtudes, and only for backward scattering angles. We present results for the 
unpolarized cross section, and for three different polarization asymmetries. We 
compare the perturbative QCD predictions for these observables with those of 
the handbag and diquark models. In order to reduce uncertainties associated 
with as and the three-quark wave function normalization, we have normalized 
the Compton cross section using the proton elastic form factor. The theoret- 
ical predictions for this ratio are about an order of magnitude below existing 
experimental data. 
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1 Introduction 



Exclusive real Compton scattering on the proton, 7p 'jp, is a promising arena for study- 
ing the short- distance structure of the proton. In the limit of large energy s/s and fixed 
scattering angle 6 in the center-of-mass frame, the real Compton amplitude should factorize 
as the convolution of a perturbative hard scattering matrix element with a nonperturbative 
light-cone distribution amplitude The distribution amplitude is for the three valence 
quarks in the proton; it describes how their longitudinal momentum is partitioned when 
their transverse separation is very small. Contributions of Fock space states with more par- 
tons in the proton's light-cone wave function should be suppressed by additional powers of 
s. However, the energy at which this asymptotic prediction of perturbative QCD (PQCD) 
becomes valid is not known a priori. Soft mechanisms such as the soft overlap (or handbag) 
model 1^, ^ and the diquark model ^ could be comparable to, or even dominant over, 
the PQCD mechanism at the presently accessible center-of-mass energies of 2-4 GeV. 

The PQCD prediction for 7p — > 7p contains a number of uncertainties. First, only 
the Born level has been computed; next-to-leading-order corrections are likely to be large. 
Related to this, the Born level prediction is proportional to a high power of the running 
strong coupling constant, [^^(yu)]^, and its renormalization-scale (/i) dependence leads to a 
large normalization uncertainty on the cross section. Second, the form of the proton distri- 
bution amplitude is not well understood. Several groups have produced model distribution 
amphtudes based primarily on QCD sum rule analyses H, ^ |10[ |ll|, [12|, |l3l. These distribu- 
tion amplitudes can lead to quite different predictions for the Compton helicity amphtudes. 
Most of the proposed distribution amplitudes tend to peak in a region where two of the 
three quarks carry relatively small fractions x of the proton longitudinal momentum. This 
has led to skepticism about the applicability of PQCD at accessible energies [|14|, ^ |^, 
because relatively soft sub-processes (relative to y/s) can reorient quarks with small x from 
the initial proton direction to the final proton direction. 

Despite all these caveats, it is still useful to know the PQCD predictions for 7p — > 7p, if 
nothing else as an asymptotic limit. There have already been four separate calculations at 



Born level [16 , 17, 0, Q9|. However, no two results agree with each other, even when the same 



proton distribution amplitudes are assumed. Given this discrepancy in the literature, and the 
need for consistent predictions from the PQCD mechanism, we undertook an independent 
recalculation of this process. Our results in fact differ from all previous work, although we 
find reasonable agreement with ref. [jl8| for a subset of the helicity amplitudes, and excellent 
agreement with ref. [l^ for forward scattering angles. 

Our results are timely in view of the experimental situation. For over twenty years, 
the highest energy wide-angle Compton data available have been from an experiment at 
Cornell pO| which investigated the energy range 4.6 GeV^ < s < 12.1 GeV^. These data 
appear to obey an approximate da/dt oc scaling law, as predicted by PQCD, although 
more precise data would be useful to confirm or refute this behavior. An experiment now 
underway at Jefferson Lab [^] should soon improve the errors on the unpolarized cross 
section and its 6- and s-dependence, in the same kinematic range as the Cornell experiment. 
This experiment also plans to measure a polarization asymmetry, the transfer of longitudinal 



2 



polarization from the incoming photon to the outgoing proton, for at least one angle. We 
shall discuss this asymmetry further in section |^. An upgrade of the Jefferson Lab electron 
beam to 12 GeV would allow for the very important extension of this experiment to 
higher energies. The proposed ELFE facility with a 25 GeV electron beam would also 
be a natural place to perform higher energy Compton measurements. 

The remainder of this paper is organized as follows. In section ^ we outline the calcula- 
tion. In section ^ we present results for the unpolarized cross section and for some different 
polarization asymmetries. Section |^ contains our conclusions. 



2 Calculation 

Since the general PQCD calculational framework for the Compton process has been described 
previously, e.g. in ref. we will be brief here. The leading-twist PQCD factorization of the 
hehcity amplitude J^^h' ^'^^ incoming (outgoing) photon helicity A (A') and proton helicity 
h {h') is given by 
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■^hh' = Y dxi dx2 dxs dyi dy2 dys S(l xj^ ^{^ -Y 

d,i •'^ j=l k=l 

x<p,ix)T,^'\x,h,X;y,h',\') m), 

where the vectors x = {xi,X2,X3) and y = (1/1,1/2,1/3) represent the quark longitudinal mo- 
mentum fractions; i labels the independent three-valence-quark Fock states of the proton, 
with distribution amplitudes 4>i{x); and d represents the sum over the diagrams that con- 
tribute to the hard-scattering amplitude Tj. 

The distribution amplitude represents the three- valence-quark component of the proton's 
light-cone wave function, after the latter is integrated over transverse momenta up to a 
factorization scale /i. (Moments of the distribution amplitude can also be defined via the 
matrix elements of appropriate local three-quark operators.) The distribution amplitude 
evolves logarithmically with /z, but (as was also done in refs. [|l^, ^) we shall neglect 

this evolution here. The full distribution amplitude for a positive-helicity proton is, in the 



notation of ref. 18 



In 33 
I^t) = TT^ / dxidx2dx3 S(l-Y,Xj) Y.(j)i{x) \i;x) , (2) 
8V6^o ^ ^ .^^ 

where 

|l;x) = \u^{xi)ui{x2)d^{x3)) , 

\2;x) = |M|(xi)rf|(x2)M|(x3)) , (3) 

|3;x) = \d^{xi)ui{x2)u^{x3)) . 

The normalization constant can be determined from QCD sum rules or lattice QCD. We 
choose /at = 5.2 x 10~^ GeV^ (as in refs. [18, |19|)- Fermi-Dirac statistics, isospin and spin 
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symmetry result in only one independent distribution amplitude, 0i; the other two are given 
by 

4>2{xi,X2,X3) = -0i(xi,X2,X3) - 0i(x3,X2,Xi) , (4) 
(f)s{Xi,X2,Xs) = 0i(x3,X2,Xi) . 

In addition to neglecting evolution of the distribution amplitude, we shall also take to be 
fixed. The Born- level cross section then scales as x f^. 

The hard scattering amplitude is computed for three collinear incoming and outgoing 
quarks. The color and electric charge dependence can be factored off of each diagram as 

tI'\x, h, A; y, h', A') = C^"^ zf^ f^'\x, y- h, A, h', A') , (5) 

where C^°'^ is the color factor, g is the strong coupling constant, and z}'^'* is the appropriate 
product of quark electric charges, while T^*^) is color and flavor independent. 

The helicities of the quarks in the hard scattering amplitude are conserved by the gauge 
interactions; therefore the proton helicity is conserved, and Ai^h' ~ h ^ h' . Parity 
and time-reversal invariance further reduce the number of independent helicity amplitudes 
to three, which we take to be 

M\\, M\\, and M\\. (6) 

In principle, 378 diagrams contribute to the hard scattering amplitude. However, 42 of 
them contain three-gluon vertices and have a vanishing color factor. Many others vanish for 
individual helicity conflgurations. 



We adopted the technique in ref. [0 of using the parity symmetry (denoted £ therein) 
between certain classes of diagrams to reduce the number that had to be computed, while re- 
serving the time-reversal symmetry as a check. All diagrams were computed by two indepen- 
dent computer programs, both based on the formalism outlined in ref. [Q. These expressions 
were found to be identical to those used in the two most recent computations 0, [1^ .[[Thus 



we agree completely with refs. [^, |T^ on the hard scattering amplitude Tj. 

The next step is to perform the four- dimensional integration in eq. (|l|) over the indepen- 
dent quark momentum fractions. For the various model distribution amplitudes |7|, p|, p|, p!0[1 
we used the coefficients of 0i listed in Table I of ref. |]18[ (and eq. (6) of ref. [|ll])- Many 



diagrams include denominators that vanish inside the (x, y) integration region, due to the 
presence of an internal quark and/or gluon that can go on shell. This is not a true long- 
distance singularity, and all the integrals are finite, diagram by diagram, but it is a technical 



obstacle to obtaining a reliable value for the integral. In the notation of ref. [|I^ , the Feynman 
is prescription leads to singular denominators of the form 

m6{{x,y)), (7) 



[x,y) + ie (x, y) 



■'■We compared our results for each diagram to the formulae given in Tables III and IV of ref. These 
tables contain three errors (found by M. Vanderhaeghcn |^ as ■well as us) in addition to one in diagram 
All that ■was published in an erratum. Ho^wever, all these errors are typographical and do not affect the 
numerical results in that paper |2^ . The errors are: In the denominator of T^^'^'^\x,'\ , [^y,} ^ [), {x^^xi) 
should be {x^,yi)] T'''-^''^)(x, t, i; y, Ti i) should be multiplied by 1/c; and the diagram related to C77 by 
ToE should be FlT, not Tii. 
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where P stands for principal part and {x,y) = x{l — ys^) — yc^, with s = sin(^/2) and 
c = cos(^/2). Diagrams can be classified by the number of singular factors found in the 
denominator; for the Compton process this number can be 0, 1, 2 or 3. The presence of 
on-shell partons in the Born-level hard scattering amplitude (for particular values of {x, y) ) 
leads to large phases in the PQCD amplitude [0, |T6|, |18[. This is in contrast to the handbag 
model, which predicts an imaginary part that is small and beyond the accuracy of the model. 

At least four different numerical methods have previously been applied to handle the 
singular integrations. Ref. |16|] performed a Taylor expansion of the numerators of the 



integrand symmetrically about each singularity. Ref. kept the e in eq. (|^ explicit, 
and evaluated the integrals for a sequence of e values tending to zero, looking for stable 



results. Ref. ||T8[ handled the imaginary parts of the singular integrals by solving the 5- 
function constraint explicitly, and carried out the real, principal-part integrals by folding 
the region of integration over at the singularity, so that the integrand is manifestly finite. 



Finally, ref. |]T9[ deformed the (x, y) integration contour into the complex plane, an elegant 
technique that requires relatively little bookkeeping for its implementation. 

We adopted a variation of the contour deformation technique . We first let 



xi=ii. a;2 = (l-^i)(l-6), X3 = (1-6)6, (8) 
yi=Vi, 1/2 = (1 -^/i)(l -r?2), y3 = {l-Vi)V2, 

so that the four independent variables (6, 6? ''7i) "'72) were integrated on the interval [0,1]. 
We then deformed the single variable rji into the complex plane, so that it ran either over 
the piecewise linear contour ^ ie 1 + ie 1, or over a semi-circular contour extending 
from to 1. Note that this simultaneously deforms both yi and ?/3, towards opposite sides 
of the real axis, while xi and X3 remain real. Inspection of the denominator factors in 
Tables III and IV of ref. JlSl shows that this deformation is sufficient to correctly bypass 
the singularities in every Compton diagram. For example, the denominator of diagram A16 
includes the factors [(a;i, yi) + ie][{x3, yi) + ie][{ys, X3) -|- ie], where Xi = 1 — Xi, yi = 1 — y^. 
Using the identity {x,y) = {y,x), the singular factors can be rewritten as [(1 — yi,xi) + 
ie][{l — yi,X3) + ie][{y3,X3) + ie], which shows that yi and ys should indeed be deformed in 
opposite directions. If the diagram happens to contain denominators of the form {xi,yi) or 
{y3,Xi), instead of {yi,Xi) or (xi^y^), as does diagram A16, then the imaginary part should 
be multiplied by an overall minus sign (or equivalently, the contour should be deformed in 
the opposite direction with respect to the real axis). 

After making these contour deformations, the real and imaginary parts of the complex 
integrals were performed separately using the Monte Carlo integration routine VEGAS |2^ . 



Two independent versions of the contour integration were implemented numerically, with 
two different choices of contour (piecewise linear vs. semi-circular), and we also varied the 
deformation parameter e, obtaining stable results. 

A third version of the integration program was constructed, which employed the Gauss- 



Legendre formalism p8| with ten points per integration variable, instead of VEGAS. Al- 
though the Gauss- Legendre errors were larger than the VEGAS errors, the two sets of results 
were completely consistent with each other (and were both inconsistent with results from 



previous work; see section 
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We carried out further checks on our integration routines. For diagrams with only one 
singular factor in the denominator, one can integrate the imaginary part analytically. Using 
this procedure we checked the imaginary part of all diagrams with one singularity. One can 
also check the diagrams with no singularities in the same manner. A second check employed 
the identity 

^ = ^-y' _ 1-^^' . (9) 
{x,y){x,z) c^{y-z){x,y) c^{y-z){x,z) 

Using eq. (||) one can reduce all three-singularity diagrams to two singularities. (These 
expressions can be reduced no further, though, because the four remaining singular variables 
are all different.) One can also reduce all diagrams that initially had two singularities to 
one-singularity diagrams, allowing their imaginary parts to be computed analytically. Our 
integration techniques were robust against all of these tests. 



Finally, Table V of ref. [IE gives detailed results for diagram A51, which has two denom- 
inator singularities. We agree completely with these results, for both the real and imaginary 
parts. We note that ref. |^ also attempted to evaluate this diagram by implementing the 



explicit e ^ method of ref. , but they obtained very different results for the imaginary 



part, compared with the results of their folding method. Ref. [0 claims that the explicit 
e — > method is not numerically stable. Since we agree with their results for diagram A51, 
we do not have cause to disagree with this claim. 



3 Results 

3.1 Comparison with previous work 

We computed the Compton helicity amplitudes for a variety of distribution amplitudes, 
which we refer to as CZ 0, GS |[, KS |], COZ HET |n|, and ASY (the distribution 
amplitude for asymptotically large energy scales, 0i(a;i, 0:2, xs) = 120a;iX2a;3). The CZ, 
KS and COZ distribution amplitudes, which satisfy the constraints imposed by QCD sum 
rules d^, are qualitatively similar. They feature a peak in 0i for xi ~ 1, ^2,3 ~ 0; that is, 
the u quark with the same helicity as the proton carries most of the momentum. The GS 
distribution amplitude has a peak in 0i for xi^s ~ 1/2, X2 ~ 0; thus it splits the momentum 
more equitably between the two quarks carrying the proton's helicity. The HET distribution 
amplitude is intermediate in shape between GS and the {CZ,KS,COZ} class. 

Before discussing our full results, we present a comparison of results in the literature. 
Here we choose the COZ distribution amplitude, since it was employed in four of the five 
existing calculations. (Only the earliest calculation |T^, which was later superseded [|1^], did 



not use the COZ distribution amplitude.) The overall cross-section normalizations in the 
literature are sometimes difficult to determine, due for example to unspecified choices for 
as- Therefore we choose to compare results for the following (normalization-independent) 
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Figure 1: Four different calculations of the polarization asymmetry An defined in eq. (10), for 
the COZ distribution amplitude. The dotted line is from ref. 
the dot-dash line from ref. ||l^, and the solid line from this work. 



17 1, the dashed line from ref. 



initial-state helicity correlation |19| , 



All = 



da^ da^ 

dt dt 

da^ 



da^ 



(10) 



where da^/dt is the differential cross section for a helicity h proton scattering off a helicity 
A photon. 

Figure |1| shows that none of the four calculations of All agrees completely with any other. 
The only two results that are very close are ours and that of ref. [0. These two curves are 
in excellent agreement for 6 < 110°; however, we do not reproduce the prominent dip of 
ref. [0 in the backward region. This statement is true for the four distribution amplitudes 
we have compared: KS, COZ, CZ and ASY [0, Figure 14 of the second reference in 



shows that the dip in All derives from ^(7tPT 7P) ^ l-^|tP + lA^tlP, and not from 



da 



dt (7iPT ~^ IP) ^ |A^|jP + Indeed, we agree with their 7|P| — >• 7p cross section for 

all angles to better than 10%, up to an overall normalization factor which can be accounted 
for by different choices for a^. We agree with the 7|j9| — ^ 7p cross section only for 6 < 110°, 
however. This suggests that the discrepancy with ref. 



hehcity amplitude A^jj. 



191 is predominantly from the single 
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The curve from ref. |]T^ has the same general shape as ours, but is offset from it. The 
phases of the dominant hehcity amphtudes given in ref. [|l^ actually agree quite well with our 
results in figs. below; the magnitudes are offset by relatively angle-independent factors. 



Ref. [18 1 finds a very large asymmetry. We have made a detailed comparison of our COZ 



results with those of ref. |]T8[, for the real and imaginary parts of the three independent 
helicity amplitudes. Each amplitude has been further split into four pieces ||2^ , according to 
the number of singular propagators in the diagram (as determined from Tables III and IV 
of ref. |jl8[ ) . The zero propagator terms (which were integrated analytically by both groups) 
agree to high precision (6 digits). The one propagator terms agree to within VEGAS errors, 

hv 



except for the imaginary part of one helicity amplitude (A^jj) which is within 10%. For 



the two propagator terms, we are in agreement on the real part of A^jj and A^||, but have 
a large discrepancy in the imaginary part. Strangely enough, for we agree on the 

imaginary part but disagree on the real part! For the three propagator terms, both the 
real and imaginary parts disagree for all three helicity amplitudes. The bulk of our overall 
numerical disagreement comes from the two propagator terms contributing to the imaginary 
part of A^jj. The two propagator terms are often 100 times larger than ours, and they drive 
ref. [0's values for ImTWiI to be roughly a factor of 10 larger than ours. 



TT 

We also calculated All for the COZ distribution amplitude using Gauss-Legendre inte- 
gration instead of VEGAS. The result agrees with our VEGAS result shown in fig. ^ (albeit 
with larger errors), and it disagrees with the other results, in particular that of ref. [|l^ for 
e > 110°. 



3.2 Helicity amplitudes and unpolarized cross section 

In figs. |]-^ we display our results for the polarized differential cross sections, 

s'^ = ^\M^.U\ (11) 

at lOTT 

for the three independent helicity configurations. Each figure plots the results for five dif- 
ferent distribution amplitudes. (For HET we shall only plot the unpolarized cross section.) 
These plots were made for a~j^ = 137.036, as = 0.3 and Jn = 5.2 x 10~^ GeV^, so they 
can be compared directly with ref. [jl8l. The phases of the helicity amplitudes are plotted 
in figs. ^|-|^; the GS distribution amplitude has a much different behavior and is therefore 
plotted separately, in fig. The phases are generally large; indeed is almost pure 

imaginary (except for the GS distribution amplitude). For reference, we also provide in Ta- 
ble 1 our numerical results for the real and imaginary part of A^jj, for the COZ distribution 
amplitude, including errors from the VEGAS integration. 

Figure |^ shows our predictions for the unpolarized differential Compton cross section, 
given by 

s^^ = ] E (12) 

^ \,\',h,h' 



along with the experimental data from ref. |2y]. For the values used as = 0.3, /n = 
5.2 X 10~^ GeV^, the predictions lie at least an order of magnitude below the data. Since 
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B (deg) 



Figure 2: The cross section for 7|P| 7|P| for five different distribution amplitudes, CZ, COZ, KS, 
GS and ASY. The results for the asymptotic distribution amplitude (ASY) have been multiplied 
by 100. 



9 



10^ 



10" 



10 



-2 _ 



10" 



/ /COZ 




ASY_ 



100 
d (deg) 



Figure 3: The cross section for ^^p^^ — ^ip^. 
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Figure 4: The cross section for 7j^p| — ^ip^. 
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Figure 6: Phase of the hehcity amphtude for 7-|-p-|- — 7jP|. 
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Figure 7: Phase of the hehcity amphtude for 7j^p-|- — 7jP|. 
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Figure 8: Phase of the three independent helicity amplitudes for the GS distribution ampHtude. 
The arrows correspond to the photon heUcities A, A' in the amphtudes A^^^ . 
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14 



d (deg) 


10^ s2Re(A<{|) 


10^ s2 Im(A^[{) 


20 


-74920 ± 240 


29200 ± 230 


30 


-15720 ± 110 


5133 ± 46 


40 


-5255 ± 15 


1301 ± 14 


50 


-2371.2 ± 8.0 


348.6 ± 5.8 


60 


-1273.6 ± 4.3 


42.2 ± 3.5 


70 


-768.8 ± 2.3 


-72.1 ± 2.3 


80 


-511.2 ± 3.1 


-115.4 ± 1.4 


90 


-369.8 ± 1.2 


-139.0 ± 1.0 


100 


-278.3 ± 1.0 


-152.03 ± 0.91 


110 


-222.4 ± 1.2 


-165.53 lb 0.90 


120 


-179.54 ± 0.70 


-183.15 ± 0.95 


130 


-144.2 ± 1.0 


-211.6 ± 1.0 


140 


-107.80 ± 0.91 


-257.1 ± 1.3 


150 


-52.9 ± 2.7 


-324.9 ± 2.6 


160 


75.9 ± 3.1 


-415.5 ± 5.0 



Table 1: The real and imaginary parts of the helicity amplitude A^jj for the COZ distribution am- 
plitude (multiplied by in units of GeV"^). The errors are from the VEGAS numerical integration. 
The values used for /jy, aem; and Ug are the same as in the rest of the paper. The normalization 
is the same as in Table V of ref. ]18U (which we found quite useful) . 



the PQCD cross section scales like , accommodating a factor of 10 by changing would 
require ~ 0.5. While this is not out of the question, and while some variation in /^y 
could be considered as well, this may be pushing the validity of perturbation theory. On 
the other hand, the shape of the curves (i.e., ignoring the overall normalization) matches the 
data quite well for the KS, COZ, CZ, and HEX distribution amplitudes. 



3.3 Normalization by Ff (Q^) 

As mentioned in the introduction, the a^{^i) scaling of the proton Compton cross section at 
Born level introduces a large normalization uncertainty into the PQCD prediction. Uncer- 
tainty in /jv also contributes. Both of these uncertainties can be removed at Born level by 



considering the dimensionless ratio |2^ 



6 da^p 



2^l2 



(13) 



where FflQ"^) is the elastic Dirac form factor for the proton at space-like momentum transfer 
Q. One might also imagine normalizing the Compton cross section by the time-like proton 
form factor. At leading order in a^, the PQCD predictions in the space- like and time-like 
regions are identical RI]; however, experimentally the time- like form factor is larger by a 



factor of about two [30, 31]. Higher order PQCD corrections can in principle account for 



15 



this factor, as Sudakov effects are different in the two regions [^. The Compton scattering 
kinematics are much closer to those of the space-hke proton form factor than the time-hke 
one, at least as far as the proton is concerned. Therefore Sudakov and related higher-order 
effects are best cancelled by normalizing with the space-like form factor. 

At leading twist, Ff (Q^) is predicted to be the same as the magnetic form factor G^^j{Q'^). 



Experimentally, these are close but not identical |31|. To normalize the experimental Comp- 



ton points, we use the experimental form factor values, 
Q'Fr{Q')^Q'Gl,{Q') 



1.0 GeV" 



^ 7-15 GeV^ 



(14) 



which are representative of the region where both scaled form factors flatten out, and are 
also similar to the highest experimental values of s available in Compton scattering.| To 
normalize the theoretical Compton curves, we recalculated the proton form factor at leading 
order in PQCD, obtaining 



216 



(15) 



where 



' 2.500 X 10^ 


(CZ), 


2.505 X 10^ 


(GS), 


3.653 X 10^ 


(KS), 


2.897 X 10^ 


(COZ) 


3.303 X 10^ 


(HET) 


.0 


(ASY) 



These results, using the wave function 



(16) 



, are precisely 



which is equivalent to that in ref. 
a factor of two smaller than several previous calculations using the same wave functions |3^ 
We do not understand the origin of this discrepancy. We do agree with the normalization 
of the hard scattering amplitude and the form factor in ref. [0 (which uses, however, a 
different representation of the proton wave function than eq. (Q)). 

Figure |10] shows the Compton cross section, normalized according to eq. (0), for both 
PQCD and the experimental data. We omit the ASY distribution amplitude, since the lead- 
ing order ASY form factor vanishes. Compared with the conventionally normalized curves 
in fig. 1^, the spread between the predictions of the three qualitatively similar distribution 
amplitudes, KS, COZ and CZ, has become much smaller. The theoretical curves also lie a 
factor of 2 to 5 closer to the data. However, they still fall about an order of magnitude below 
the data at the widest scattering angles. (The HET distribution amplitude does slightly bet- 
ter than this.) Thus it seems unlikely that the elastic proton form factor and the Compton 
scattering amplitude are both described by PQCD at presently accessible energies, unless 
there are large higher-order and pro cess- dependent corrections. 



^If one equates the four- momentum transfer to the proton in the two processes — in the form factor 
and —t in Compton scattering — then the corresponding Compton s = 2(5^/(1 — cos 6') should actually 
be considerably bigger than . At 90°, for example, s = 2Q^ . Unfortunately, there are no experimental 
Compton data with s this large (all have —t < 5.3 GeV^), so there is not a good overlap with the region (|lT 
where the elastic form factor is beginning to scale properly. 
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Figure 10: The scaled unpolarized Compton cross section, normalized by the scaled elastic proton 
form factor, as in eq. (|l3|), for five distribution amplitudes, compared with the experimental data [pO| , 
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Figure 11: The initial state helicity correlation All in perturbative QCD for five distribution 
amplitudes. Also plotted is the handbag model prediction for = 4 GeV (GRV) and a 
diquark model prediction 0]. 



3.4 Asymmetries 

Various polarization asymmetries can be constructed from the helicity amplitudes. These 
observables may provide additional diagnostic power for uncovering the Compton scattering 
mechanism, beyond what the unpolarized cross section provides. 

Figure |Tl| presents the perturbative QCD results for the initial state helicity correlation 
All defined in eq. (|l3). Also shown is the handbag model prediction for = 4 GeV, 
where the form factors Rv,a were evaluated using the parton distribution functions of 
GRV |35]. In leading-twist PQCD, the proton helicity is conserved. The handbag model 



does not inherently require proton helicity conservation, but it has been assumed in ref. 



Thus the PQCD and handbag curves for All in fig- O can be equated to the longitudinal 
photon-to-proton polarization transfer asymmetry, which is slated to be measured for at 
least one scattering angle in an upcoming experiment [^. The diquark model analyzed in 
ref. 1^, P] has nonvanishing proton helicity-flip amplitudes at finite s, making All and the 
polarization transfer into distinct asymmetries. We plot the diquark prediction for All from 
ref. Figure |ll] shows that PQCD gives quite different qualitative behavior from both 
the handbag and diquark models for All, and they should be distinguishable with the help 
of experimental data at just a couple of backward scattering angles. A caveat is that the 
GS curve is somewhat oscillatory, so one might wonder whether a distribution amplitude 
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Figure 12: The photon spin transfer coefficient Dll in perturbative QCD for five distribution 
amphtudes. Also plotted is the diquark model prediction for = 4 GeV (standard DA) Q. The 
handbag model predicts Dll = 1- 



'between' GS and the {CZ,COZ,KS} class of amplitudes could produce behavior similar to 
the handbag model. 

One can also define a photon spin transfer coefficient 



dt dt 



where now da^'^' /dt is the differential cross section for initial and final state photon helicities 
A and A', and unpolarized incoming and outgoing protons. Figure |T2| gives the PQCD 
predictions for this asymmetry, as well as that of the diquark model for = A GeV and 
a 'standard' distribution amplitude 0. The handbag model predicts Dll = 1, basically 
because the helicity-fiip quark Compton amplitude 7|g — > 7j^g vanishes at Born level for 
massless quarks. 

The final asymmetry we plot is the photon asymmetry 

da± _ da\\ 
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Figure 13: The photon asymmetry T, in perturbative QCD for five distribution ampHtudes. Also 
plotted is the diquark model prediction for = A GeV (standard DA) The handbag model 
predicts S = 0. 



20 



where da±/dt and da\\/dt are the differential cross sections for hnearly polarized photons, 
with the polarization plane perpendicular or parallel (respectively) to the scattering plane. 
Generation of this asymmetry requires a nonzero photon helicity-fiip amplitude; hence the 
asymmetry vanishes in the handbag model. Figure |13| plots the PQCD and diquark pre- 
dictions. The diquark prediction is shown for = 4: GeV and a 'standard' distribution 
amplitude; for another distribution amplitude E can become positive in the backward region 
instead of negative H. This asymmetry has actually been measured |PB[, however only for 



= 3.45 GeV and cos 9 > 0.8. A high-energy wide-angle measurement would be very 
useful for distinguishing between handbag and PQCD mechanisms. 



4 Conclusions 

Motivated by conflicting results in the literature, we have recalculated the fixed-order. Born 
level predictions of perturbative QCD for proton Compton scattering, for five different dis- 
tribution amplitudes. While our results do not agree with those of any previous group, they 
do agree very well with those of ref. [jl9[ for 6 < 110°, and the differences for 6 > 110° seem 



to be dominated by a single helicity amplitude, Al||. 

From the helicity amplitudes we computed three separate polarization asymmetries. Ex- 
perimental measurements of these asymmetries could be used in conjunction with the un- 
polarized differential cross section in order to help shed light on the mechanism involved in 
the Compton scattering process. 

We also have attempted to reduce the uncertainty in the overall normalization of the 
Compton cross section by normalizing it by the square of the elastic proton form factor. 
This exercise reduces the spread in the theoretical predictions, but it leaves them an order 
of magnitude below the data. Unfortunately, this result makes it difficult to simultaneously 
explain the current data on the elastic proton form factor and on Compton scattering in terms 
of perturbative QCD, without appealing to large uncalculated higher-order and process- 
dependent corrections. 

Acknowledgments 

We thank Stan Brodsky for suggesting this work, and for useful discussions and com- 
ments on the manuscript. We are grateful to Marc Vanderhaeghen and particularly Andreas 
Kronfeld for detailed discussions of their work and for providing us with unpublished data 
files. We thank Markus Diehl for several helpful remarks and suggestions. 

References 

[1] G.P. Lepage and S.J. Brodsky, Phys. Rev. D22, 2157 (1980). 

[2] A.V. Radyushkin, Phys. Rev. D58, 114008 (1998) [|hep-ph/9803316 |. 



[3] M. Diehl, T. Feldmann, R. Jakob and P. KroU, Eur. Phys. J. C8, 409 (1999) jEi^ 
ph/98112"53| . 



21 



[4] M. Diehl, T. Feldmann, R. Jakob and P. Kroll, Phys. Lett. B460, 204 (1999) jEi^ 



ph/9903268 



[5] 
[6] 
[7 



[9 

[lo; 
[11 

[12 

[13: 

[14 

[is; 

[16 
[IT 

[is; 

[19 



[2o; 

[21 



p. Kroll, M. Schurmann and W. Schweiger, Int. J. Mod. Phys. A6, 4107 (1991). 

C.F. Berger, Diplomarbeit, Technical University of Graz (November 1997). 

V.L. Chernyak and I.R. Zhitnitsky, Nucl. Phys. B246, 52 (1984). 

M. Gari and N.G. Stefanis, Phys. Lett. B175, 462 (1986); 
Phys. Rev. D35, 1074 (1987). 

I.D. King and C.T. Sachrajda, Nucl. Phys. B279, 785 (1987). 

V.L. Chernyak, A.A. Oglobin and I.R. Zhitnitsky, Yad. Fiz. 48, 841 (1988) [Sov. J. 
Nucl. Phys. 48, 536 (1988)]; Z. Phys. C42, 569 (1989). 

N.G. Stefanis and M. Bergmann, Phys. Rev. D47, R3685 (1993) ||hep-ph/9211250|] . 

M. Bergmann and N.G. Stefanis, Phys. Lett. B325, 183 (1994) ||hep-ph/9403^ |. 

R. Eckardt, J. Hansper and M.F. Gari, Phys. Rev. D50, 26 (1994); 
|hep-ph/9607380| . 

N. Isgur and C.H. Llewellyn Smith, Phys. Lett. B217, 535 (1989). 
J. Bolz and P. Kroll, Z. Phys. A356, 327 (1996) [|hep-ph/9603289| . 
E. Maina and G.R. Farrar, Phys. Lett. B206, 120 (1988). 

G.R. Farrar and H. Zhang, Phys. Rev. D41, 3348 (1990), D42, 2413(E) (1990); Phys. 
Rev. Lett. 65, 1721 (1990). 

A.S. Kronfeld and B. Nizic, Phys. Rev. D44, 3445 (1991), D46, 2272(E) (1992). 

M. Vanderhaeghen, P.A. Guichon and J. Van de Wiele, DAPNIA-SPHN-96-36, Pre- 
sented at Workshop on Virtual Compton Scattering, Clermont-Ferrand, France, 26-29 
Jun 1996; Nucl. Phys. A622, 144c (1997); 

P.A. Guichon and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 41, 125 (1998) fhep- 
ph/9806305|| . 



M.A. Shupe et ai, Phys. Rev. D19, 1921 (1979). 
A.M. Nathan, |hep-ph/9908522| ; 

J.P. Chen et al. [Hall A Collaboration], PCCF-RI-99-17. 



[22] See e.g. [http : //www, cebaf . gov/div_dept /phys ics_division/GeV. html 



22 



[23] ELFE "blue book", J.-M. DeConto, ed., Institute des Sciences Nucleaires de Grenoble, 
RME02 (October 1993); 

R. Brinkmann et al, Nucl. Phys. A622, 187c (1997); 

R Hoyer, Nucl. Phys. A622, 284C (1997) ||hep-ph/9703"46^ ; |hep-ph/9911486i 
G.R. Farrar and F. Neri, Phys. Lett. 130B, 109 (1983), 152B, 443(E) (1985). 
M. Vanderhaeghen, private communication. 
A.S. Kronfeld, private communication. 
G.R Lepage, J. Comput. Phys. 27, 192 (1978). 

W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical Recipes in 
C: The art of scientific computing (Cambridge Univ. Press, 1996). 

S.J. Brodsky, private communication. 

M. Ambrogiani et al. [E835 Collaboration], Phys. Rev. D60, 032002 (1999). 
A. F. Sill et al, Phys. Rev. D48, 29 (1993). 

T. Gousset and B. Pire, Phys. Rev. D51, 15 (1995) [|hep-ph/9403293 |. 

V. L. Chernyak, A. A. Ogloblin and I. R. Zhitnitsky, Z. Phys. C42, 583 (1989); 
C. Ji, A. F. Sill and R. M. Lombard, Phys. Rev. D36, 165 (1987); 
N. G. Stefanis, Phys. Rev. D40, 2305 (1989). 

C. E. Carlson and F. Gross, Phys. Rev. D36, 2060 (1987). 

M. Cluck, E. Reya and A. Vogt, Z. Phys. C67, 433 (1995); 

M. Cluck, E. Reya, M. Stratmann and W. Vogelsang, Phys. Rev. D53, 4775 (1996) 
|hep-ph/ 950834711 . 



[36] G. Buschhorn et al, Phys. Lett. 37B, 211 (1971). 



23 



